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OO . Abstract 

in ■ 

£NJ ■ We present the results for two-loop MSSM corrections to the relation between pole and 

running masses of heavy quarks up to the order 0{o? s ). The running masses are defined in the 
\ DR scheme, usually used in multiloop calculations in supersymmetric theories. The analysis 

of the value of these corrections in different regions of the parameter space is provided. 

C^' 1 Introduction 

<d : 

^Pi \ Recently a lot of theoretical and experimental efforts have been spent to find some new physics 
.£h ! beyond that described by the Standard Model (SM). One of the most popular ways to extend SM is 
^ ■ to make our world supersymmetric at some higher scale than those currently tested by experimental 
facilities over the world. A way of doing it minimally and at the same time to have renormalizable 
local quantum filed theory is to consider the N = 1 SUSY field theory. It is obtained from SM 
by replacing SM fields with groups of fields, forming representations of the N = 1 SUSY algebra, 
and introducing an additional Higgs doublet required by supersymmetry. This way we obtain the 
Minimal Supersymmetric Standard Model (MSSM) [DGOIIII 

Huge theoretical research in this direction was made in the last few years as well as parameter 
space of this model was heavily constrained by experimentalists. Today people consider radiative 
SUSY corrections to SM precision observables and try to observe SUSY particles directly in an 
experiment. The precision with which masses of heavy quarks are known at present, requires 
that higher order corrections, in addition to the leading one-loop SUSY corrections, be taken into 
account. 

It is the aim of the present article to calculate the "leading" 0(a 2 s ) MSSM corrections to the 
pole masses of heavy quarks and to see how they affect the SUSY particle spectra obtained from the 
renormalization group analysis with universal boundary conditions at the GUT scale. In fact we will 
be interested in the relation between the pole mass of a heavy quark and the running heavy quark 
DR mass. The latter is defined as a quantity computed in dimensional reduction and renormalized 
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minimally. 1 By "leading" corrections we mean those, which are mediated by strong interactions 
and formally they should be dominant. However, as it is known from one- loop calculations, in the 
case of 6-quark there exists a large contribution from a stop-chargino loop enhanced by the t-quark 
Yukawa coupling, tan/3 or by the supersymmetric Higgs mass parameter ji. So, in this sense, our 
results for the 6-quark are not complete. But, they can be used to see the typical size of these 
corrections and already make some conclusions on whether they should be accounted for in an 
accurate analysis or not. 

To compute needed diagrams we use the large mass expansion procedure and for the moment we 
restrict ourselves only to terms up to the C(m^/M| USY ) order in the case of 6-quark and up to the 
(9(m^/M| USY ) order in the case of t-quark. It is a good approximation for the 6-quark. However, 
as our numerical analysis shows, such approximation does not work well in the case of t-quark and 
higher terms in large mass expansion should be taken into account. We suppose to calculate missing 
corrections from stop-chargino loops to the 6-quark pole mass and to provide more terms in the 
expansion in a relation between the top pole and DR masses elsewhere. 

The paper is organized as follows. In section 2, we define quantities we want to compute and 
make comments on the choice of our renormalization scheme. In section 3, we describe our model 
- the supersymmetric QCD (a subset of MSSM, relevant to the calculation of a™ corrections) 
and present in detail the renormalization procedure. Section 4 contains our results and numerical 
analysis of the effect produced by our correction on supersymmetric particle spectra. The latter 
were computed with the help of the SoftSUSY program [S| under universal boundary conditions at 
GUT scale. Finally, section 5 contains our conclusion. 

2 Pole mass and choice of renormalization scheme 

First, we would like to note that there are several quark mass definitions. This is mainly due to 
the fact that quarks were not ever observed as free particles. Therefore, the definition of their 
masses rely heavily on theoretical constructions. Different definitions exist referring to different 
renormalization schemes used in quantum field theories. We will mention only a few most popular 
ones — MS mass m, pole mass M and DR mass m. In the nonsupersymmetric QCD the relations 
between these three masses are known (to two-loops between DR and pole masses and to two 
[HI Ej and to three [THl HI] loops between MS and pole masses) and in the supersymmetric QCD case 
there exists only one- loop relation between DR and pole masses of heavy quarks [T2*| IT3*| IT4*| HH| ITH] . 
In this work, we will establish the two-loop relation between DR and pole masses in supersymmetric 
QCD. 

In DR scheme, particle masses and couplings depend explicitly on the renormalization scale jx 
which is taken in applications equal to the characteristic scale of a studied process. The renor- 
malization scale dependence of these quantities could be conveniently described by renormalization 
group equations. They are known very well for MSSM and supersymmetric QCD. It should be noted 
that DR mass is a short distance quantity. That is, it is sensitive only to short distance effects. To 
describe processes with the characteristic scale of the order of quark mass itself, a different on-shell 
mass definition is used and here pole mass of a particle comes into play. 

The particle pole mass is defined by the singularity of the corresponding two-point function. As 
explicit perturbative calculations showed, the pole mass of a quark is an infrared finite and gauge 

1 The DR scheme which is similar to the MS scheme, is usually used for calculation in supersymmetric models, 
because contrary to the MS scheme it manifestly preserves supersymmetry, at least for not very high loop orders 5 . 
The details about the renormalization scheme used will be given below. 
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invariant quantity. For this reason it is considered as a physically meaningful quantity fZj . We will 
restrict ourselves to perturbation theory only and will not analyze the exact nature of the two-point 
function singularity which may involve some nonperturbative dynamics. In the present paper, we 
employ the definition of the pole quark mass, where it equals the value of p, at which the inverse 
quark propagator turns to zero. 

The pole mass M of a quark is defined as a formal solution for p (in the Minkowski metric) at 
which the reciprocal of the connected full propagator equals zero 

p — m — E(p, m) = 0, (1) 

where E(p, m) = m'Ei(p 2 ,'m) + (p — m)S 2 (p 2 ,m) is the one-particle-irreducible two-point function, 
m may stand for the bare or renormalized mass, m-bare o r m ren , depending on the prescription used 
in evaluating S. The solution to eq. (JTJ) is sought order by order in perturbation theory. To two 
loops we have 

M = m + S (1) (m,m) + S (2) (m,m) + S (1) (m,m) S (1) '(m,m) + C(S (3) ), (2) 

where is the L-loop contribution to X, and the prime denotes the derivative with respect to 
the first argument. In what follows, we will be interested in the relation between the pole quark 
mass M and the running DR mass m computed in MSSM up to the 0(a 2 ) order. Technically, to 
solve this problem, we need to evaluate ~ 200 two-loop propagator type diagrams involving many 
different mass scales (see Fig. [Q. In general, it is quite a complex problem to be solved exactly. 
However, in our case all mass scales can be divided into two groups denoted by m so{t and m hard , 
such that scales from m so f t are much less than each of the scales from the mh ar d group 

case of 6-quark: m& m t , rriq, mg 
case of t-quark: m t rriq, mg 

Here mg is the gluino mass, rriq stands for different squark masses. We also explicitly denoted which 
mass scales belong to soft or hard groups in the cases of b- and t-quarks. 

Given this mass hierarchy, one can employ the large mass expansion procedure to reduce evalu- 
ation of multi scale two-loop integrals to the calculation of single scale on-shell two-loop integrals, 
two-loop tadpole integrals with two scales and products of one-loop on-shell integrals and one-loop 
tadpole diagrams. To compute two-loop single scale on-shell integrals, we made use of the ON- 
SHELLY package ^H] • For evaluation of two- loop tadpole diagrams with two scales, the recurrence 
relations of 19, were used. Calculation of one- loop self-energies, including their derivatives with 
respect to momentum, is an easy task and we will not make further comments on it. 

Now let us make some comments on the choice of a renormalization scheme. Here we use the same 
renormalization prescription, as in [7j. To be more specific, we use regularization by dimensional 
reduction, a modification of the conventional dimensional regularization, originally proposed in j^j. 
In this procedure, the vector and spinor algebras in the numerator of Feynman diagrams are four- 
dimensional, which is the requirement imposed by supersymmetry. However, to make sense from 
divergent momentum integrals we need to regularize them by noninteger space-time dimension 
d = 4 — 2e. For quantum corrections not to break gauge invariance, we need a cancellation of 
squared momenta in both the numerator and the denominator of Feynman integrals. Thus, the 
momenta should form d- dimensional subspace in four dimensions. As the momenta become d- 
dimensional, four-vectors naturally split into true d-vectors and so-called e-scalars, obtained in 
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the process of dimensional reduction of the original four- dimensional Lagrangian. The e-scalars are 
nothing, but matter fields. Their appearance is the only difference from the conventional dimensional 
regular ization. 

We would like to note that, in general, renormalization of e-scalars and their interactions is not 
identical to that of vectors, so the original four-covariance may be spoiled by quantum corrections. 
Moreover, quantum corrections may also generate a mass for e-scalars. There is arbitrariness in 
choosing the renormalization scheme for this mass (20] • A consistent way is to choose the finite 
^-scalar mass counterterm, so that the pole (and renormalized) mass of the e-scalars equals zero. 
The e-scalar field renormalization is left minimal. Other renormalizations are done minimally, that 
is by subtracting only poles in e. 



3 Definition of the model and renormalization procedure 

For our purposes here we do not need the complete MSSM, but only its part — supersymmetric 
extension of QCD (SUSY QCD). In the superfield formulation SUSY QCD Lagrangian consists of 
two parts: the rigid supersymmetric part and the part containing soft supersymmetry breaking 
terms. 
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where Q = Q L , U = U R , D = D R , T a = X a /2 (a=l,...,8) with A a being the Gell-Mann matrices 2 
and the gauge field strength tensors are 



W a = -—D 2 e- 2 ^ Taya D a e 2 ^ Tav \ 
893 



W h = -—D 2 e 2 ^ Tava D^e- 2 ^ Taya 
893 



The superpotential of the ordinary SUSY QCD contains only quadratic mass terms for chiral 
(antichiral) superfields. Considering SUSY QCD as a part of MSSM, where all masses of particles 
are generated through the Higgs mechanism, the superpotential takes the following form 



y d DH{Qi + yJJHZQ* + fxH{H l 2 ] , (4) 

where the and y u are the Yukawa coupling constants carrying the generation indices which have 
been suppressed (as well as group indices). In this work we will need the Higgs fields only to 
generate masses of our particles and will not be interested in the detailed structure of the MSSM 
Higgs sector. 

To perform the supersymmetry breaking that satisfies the requirement of "softness", we intro- 
duce a gluino mass term, soft masses of scalar superpartners of quarks and soft trilinear couplings 
with Higgs fields, which may be written in terms of N — 1 superfields, provided one introduces 
external spurion superfields |2L 
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y d A d DH[Q j + yJ^UHlQ*] 6 2 + h.c. 



2n2 



9 2 9 



(5) 



2 Note that U and D are colour antitriplets and the anticolour generator is — (T a )" 
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The chiral superfields have the following component fields representations 

q = q L - t ea^q L -hem^d^ L + V26q L -^=eeea^q L + eeF q , (6) 
u = u R -iea^u R -]eeeed^d^u R + V2eu R -^ r e9ea l "d ti u R + eeF u , (7) 

4 y2 

D = d R -i9a"9d„d R - ]eem»dJ R + V26d R - -^=999a t *d ti d R + 99F d (8) 

4 \/2 

and for the vector superfields in the Wess-Zumino gauge we have 

V a = 9a^9G a ^ + 999l a + 999g a + ^9999D a . (9) 

Then, rewriting everything in the component fields, the Lagrangian of SUSY QCD takes the form 
(Qr = (uR,d R )) 



1 

— i 
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C = --G%G a ^ + il a ^D^ a + nR^D^q R + iq L G»D^ L 



+ (D»q R y D,~q R + (Z^g L ) T D^~q L - g 3 q R T a D a q R + g 3 q L T a D a q L 
+ V2g 3 q R T a ~g a q R - V2g 3 q L T a £ a q L + V2g^ R T a g a q R - V2g^ L T a g a q L 



l D a D a + l _w + l _ D , D > + g 2 ^ D % + g^h~D% + gJi'&D'h] 



2 * 2 

,0 r . , ..TP lO 



+ /^°^ 2 - Vdd R F h od L - y u u R F h ou L 



yd(F dR tiid L + d R h\F dL ) - y u (F UR h° 2 u L + u R h^F, 



"L , 



+ Vddiih'ldL + y u u R h° 2 u L + y d A d d R h^d L + y u A u u R h° 2 u L + h.c. 

- -yg g - -yg g - A%z?£ - M^urUr - Mp R d R , 
where the gauge covariant derivative and field strength are defined as 

D, = d, + ig3 T a G«, 

ria Q r~ia a r~ia rabcsib r~ic 

[iv — H v ~~ °^ix ~ 93J ■ 

Here we also added auxiliary SU(2) x £7(1) fields D l and D', which are necessary to describe correctly 
the mass generation mechanism, a 1 denote Pauli matrices and Y q is a hypercharge of a particle q. 
Eliminating auxiliary fields F q , F q , D a , D\ and D', and introducing four- component spinor notation 
similar to [TJ E2] , one gets 

C = -\G%G a ^ v + l -rrD^ a + iqYD^q + (D^q R y D^q R + (D^ D^q L 

+ V2g 3 qP L T a g a q R - V2g 3 qP R T a g a q L + V^g^T^P^ - V2g^ L T a g a P L q 

- g ^(q L T a q L - q R T a q R )(A L T a q L - q^ R T a q R ) - ^g a g a + C M , (10) 
C M = ~M~q L q L - Mlu R u R - M~d R d R - y d h\dd - yj^nu 



2/d (^i^l) - Vd h ihid R dR - y 2 u [h° 2 u L j h\u L - y 2 Ji\h\u R u R 
y d h°A d d R d L + y u h 2 A u u R u L - y d fih 2 d R d L - y u ^h\u R u L + h.c (11) 

where Pl = (1 — 75) /2 and Pr = (1 + 75) /2 are the chiral and antichiral projectors. 
After Higgsing and taking into account that 

m d = y d vi, m u = y u v 2 , tan^ = — , sin W = -5- — o, 

«i 3i + £2 

m 2 z = \{gl + gl){vl + vl) = iyl + v 2 2 ), 

z z sin u\y 



Y„ 



cos2/? = ^— \ e =11 + 13 
it is convenient to rewrite Cm from eq. (jllj) in the following form |2S] 



£ A f = -m g gg - m 2 - L q L q L - m 2 - R u R u R - m\ R d R d R - a q m q q L q R - a q m q q R q L (12) 



with 



m 



M| + m| cos 2/5 (J| - e 9 sin 2 9 W ) + m, 



2 



m 2 ^ = M~ + e u m| cos 2/5 sin 2 6^ + m 2 , a u = A u — fi cot /3 
= Mj + m| cos 2/? sin 2 6*vy + m 2 d , = A d — jj, tan /3 

where e q and are the electric charge and the third component of the weak isospin of a particle 
q, m q is the mass of the partner quark and fi is a Higgs mass parameter. The M q , M s , and Mj are 
soft SUSY breaking masses, A u and A d are trilinear couplings as in eq. (0). The family indices have 
been suppressed. 

Now we must diagonalize the mass matrix M 2 - of squarks to determine the physical mass eigen- 
states 

^(<^ : )^V(< 4)*. 03) 

According to eq. (fT3j) .M- is diagonalized by a unitary matrix 7^. Assuming that CP violating 
phases occur only in the CKM matrix, we choose TZ^ to be real. The weak eigenstates q L and q R 
are thus related to their mass eigenstates g x and q 2 by 

f^UWH K*=( COS ^ Sin M, (14) 
with being the squark mixing angle. The mass eigenvalues are given by 



• = \ (m 2 q L + m\ R =F ^(m^ - m 2 ^) 2 + 4a 2 m 2 ) . (15) 



l 9J,a 2 



By convention, we choose q 1 to be the lightest mass eigenstate. Notice that m qi < m qL R < m q2 . 
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Figure 1: Feynman diagrams for the O(a^) contribution to quark self-energies (last two diagrams at 
the one-loop level of the order 0(a s )). The diagrams in the first line are the pure QCD contribution. 

For the mixing angle 9g we require < 9„ < tt. Then one can write 



COS 6g 



-a q m q 



sin 6q 



m\ L - ml 



(16) 



cos 9q | > ^ if rriq L < rriq R and | cos^l 



< 75 if m Q R < m <iL- 



Moreover, 

Making transition in eq. ()10|) from and q~i to q~\ and g2 with the help of the matrix from 
eq. (J 14)) we get for quark-gluino-squark trilinear and four-squark vertices complicated expressions 
due to squark mixing which we do not present here. To produce input amplitudes of the diagrams 
needed in our calculations we made use of the Mathematica package FeynArts [21], which has SUSY 
QCD as part of MSSM model [25J. The Feynman diagrams contributing to quarks self-energies in 
the one and two loop order are shown in Fig. ^ 

Note again that doing these calculations we work in the minimal subtraction scheme MS adopted 
to supersymmetry — dimensional reduction with minimal subtraction DR . To obtain a relation 
between pole and DR quark masses, we need to know DR renormalization constants for quark, 
gluino, squark, e-scalar masses, gauge and Yukawa charges 3 as well as renormalization of squark 
mixing angle. We will describe evaluation of each of these quantities in the next subsections. 



5 By Yukawa charge we mean the Yukawa coupling of e-scalar to fermions 
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3.1 Renormalization of QCD sector 

First we have to renormalize the charge and the masses of quarks and fields. The relations between 
bare and renormalized parameters are given by 

« s ,o = fi 2e Z as a s , (17) 

m qfi = Z mq m q , (18) 

Go = Z l G /2 G, (19) 

* = Z l J 2 q, (20) 

where p, is a renormalization scale. 

The charge renormalization constant can be obtained from the renormalization constants of the 
fields and vertex with the use of the following relation: 

Z as = Z 2 Gq Z q 2 Z G X , (21) 

where Z q and Zq correspond to renormalization of quark and gluon fields and Z q Q q renormalizes 
the quark-gluon interaction vertex. 

One can also renormalize the gauge parameter 4 with the help of the above-introduced constant 

Z G 

& = Z G Z. (22) 

However, in the gauge invariant quantities (like, e.g., the pole mass) the gauge parameter drops out 
already from the bare expression and renormalization (J22)) is not needed. 
Our results for renormalization constants in the DR scheme read 

- 3C A - 2C F ) \ + (-6 + 3C A - 2C F ) -) , (23) 
Z q = i-£c F (£ + i)±, (24) 

(25) 

1-TtM- 1 OCj +^C a )-, (26) 




AqGq 

Z as = 1-^(3^-6)-. (27) 

For SU(N C ) we have Cp = (JV^ — 1)/(2N C ) and C A = N c . Here we explicitly put the number of 
quark flavors to six. 

There are several possibilities to check the correctness of our results. The coupling renormaliza- 
tion constant (|27j) can be compared with that obtained from the one-loop MSSM beta function for 
SU(3) group. We have also checked that in the case of ordinary QCD we recover renormalization 
constants given in [7] . The 2- loop DR quark mass renormalization constant ()23|) has been obtained 
from eq. (J2J) and contains poles in e remaining after all other renormalizations were done. Here we 
mean that one should perform one-loop renormalizations of strong coupling constant and particle 
masses. A subtle point is the renormalization of the Yukawa interaction of e-scalars and quarks. 

4 Gauge parameter £ is defined in such a way that £ = and £ = 1 correspond to the Landau and Feynman gauge, 
respectively. 
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In SUSY QCD discussed here, one does not need to introduce a separate coupling constant Y for 
this interaction, as supersymmetry insures that it is renormalized in the same way as a strong cou- 
pling constant 5 . However, to make our renormalization procedure as general as possible ( and also 
applicable in the case of nonsupersymmetric QCD for checks), we introduce an additional indepen- 
dent constant Y and renormalize it separately. So, here we follow the renormalization procedure 
described in j7j . Renormalizing a physical quantity like a pole mass, we ignore field renormalization 
constants which is a justified procedure, provided we are evaluating the contribution of countert- 
erms for the whole sum of diagrams. The 1/e 2 poles in Z mq can be checked with the use of the 
renormalization group technique. Indeed, let us write 



m 



a 




e 2 



{21 



where = EZi( a s/^Y Z {l ' n) . Now using the renormalization group equation 

d 



= ^ 2 ^2" m "' 



we obtain 



1 d (i) 

2 dg 



(l + ^ + Emr)' 1 * = (29) 



drrii I 2 dg 



(g is related to a s via g 2 /4:7i = a s ). The formulae presented are general and are also applicable in 
a theory with spontaneously broken symmetry [HE!- In the case of SUSY QCD, from the previous 
formula it follows that 

2ZW = (z^) 2 + 2f3^Z^ 1 \ (30) 

Substituting into this relation the expressions for Z 1 - 1 ' 1 ' = —2Cp and (3^ = — ^{3Ca — 6), we get 
the same result for the 1/e 2 poles as in eq. (|23|). 



3.2 Gluino mass renormalization 



The renormalization of the gluino mass is performed in a way similar to quark mass renormaliza- 
tion. The relevant diagrams with gluino-gluon and quark-squark loops are shown in Fig|21 Explicit 
calculations give the following results for the gluino mass mg i0 = Z m -mg, gluino wave function 

5 See the section about the e-scalar sector renormalization for details 

6 In fact too is not a renormalization group invariant in our renormalization scheme. However, to check the leading 
poles in e we can treat it as an RG invariant, as different prescriptions used to renormalize the e-scalar sector give 
different results only in the subleading order in e. This point will be discussed in detail in the section concerning the 
e-scalar sector renormalization. 
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g Q 
Figure 2: One-loop diagrams contributing to gluino mass counterterm. 



renormalization constants Zg and for gluino-gluon-gluino vertex renormalization constant ZgGq 

a 1 
Z m - = 1 + ^1(6-300-, 



a 1 
Z g = 1-^(^ + 6)-, 



(31) 
(32) 
(33) 



From eqs. (J25J), (jS2j) and (|3l?|) one can extract the gauge coupling renormalization constant 

Za s = Z\ G ~Z~ 2 Z G 1 , (34) 



which, of course, coincides with eq. fl27|) . This serves as an additional check. The result for the 
gluino mass renormalization constant can also be verified with help of the known gluino mass beta- 
function. 



3.3 Squark sector renormalization 

To discuss the renormalization of the squark sector let us start with the Lagrangian written in terms 
of physical states. The bare Lagrangian in terms of bare parameters and bare fields is given by 

4^ = d^d^-ml^ + C^, (35) 
£j g f = -V2g 3 q(P R cose q - P L sinB q )T a ~g a q x 

+ V2g 3 q (P R sin 6~ q + P L cos 9 q ) T a g a q 2 

- V2g 3 ^T a g a (P L cos Og-Pn sin O^q 

+ V2g^ 2 T a g a (P L sin 6 q + P R cos 6 q )q (36) 
We can rewrite this Lagrangian in terms of renormalized physical quantities and counterterms 

^ = C q + 5C q (37) 

5Cg = dZ^d^d^qi - m-.^) - Jm-.g^ - bm\ 2 (q 2 q x + q x q 2 ) + 5Z m C qM (38) 

where C q and C q g q have the same form as (|33jl and (|3^jl . It should be noticed that this set of 
counterterms is sufficient to cancel all rising divergencies. The 5fh\ 2 counterterm is used to render 
nondiagonal squark self-energies (i ^ j) finite. In practical calculations, it is convenient 

7 We do not include the bare Lagrangian for a four-squark interaction, because the counterterms associated with 
the four-squark vertices will give contribution of the 0{a' i s ) order. 
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't \ J Q '/ \o J Q 



Q g 
Figure 3: One- loop diagrams contributing to the squark mass and squark mixing angle counterterm. 

to trade 5m\ 2 for the counterterm for the squark mixing angle. Indeed, we can diagonalize the 
Lagrangian (}3*Tj) with the following field redefinitions: 



<?! \ _ / 1 56~ q \ / qi 
q'J~\ ~56 q 1 )[q 2 



(39) 



where we have taken into account that 56 q ~ 0(a s ), cos (56 q) = \+0(a 2 s ) and sin (56 q) = 56q+0(a 
Demanding that the new mass matrix of squarks 

1 56- q \ / m\ + 5Z~ q m\ + 5m\ 5m\ 2 \ / 1 -56 q 

-56 q 1 J V 5m\ 2 ml + 5Z^n% + 5m\ ) \ 56 q 1 

must be diagonal and expanding eq. (j4Ti|) in the gauge coupling constant a s , one gets 



3^ 



(40) 



56 q = /™' 2 2 = , (41) 

m\ x - m| 2 m\ x - m| 2 

where the subscript "div" stands for the 1/e part of the expression. Rewriting eq. (|57j) in terms of 
the new fields g-, we come to the following renormalization prescription for the squark sector: 

m lo = m l + Sm l = z m*™\, ( 42 ) 
%o = (l + \Zq)q i = Z 1 J 2 q i (43) 

%0 = 6q + 56q = Z 6 ,6q (44) 

Explicit calculation of squark self-energies gives (see diagrams in Fig. |3J) the following expression 
for the renormalization of the squark mixing angle 



sin (4:6 q ) 4m 9 m g cos (26 q ) \ 1 



56 q = C F ^ + ^Zirip _. (45) 

4tt V 2 - m| 2 / £ 

Note that renormalization (J45|) of the mixing angle is minimal. It is the correct choice, if one wants 
to have the final result expressed in terms of the running DR sparticle masses. 

Summarizing the result of calculation, the renormalization constants for the squark sector read 

a 1 

Z q = 1-C^(e-1) ? (46) 



= 1 + ^y- ^ ^ (m~ 2 - m| ) sin 2 26 q + 4 (m ? mg sin 26 q - m\ - m^j ^ , 



(47) 
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Figure 4: One-loop diagrams contributing to the e-scalar wave function and mass counterterms 



Z m ^ = 1 + ^ ( ( m l ~ m l ) sin2 26 q - 4 (m q m g sin 26 9 + m\ + ml) j ~ , (48) 

Zm = + ( 49 ) 

^ = i + f 1 - (50) 

ca- 
using eqs. (JS21) (0E|) and (jH)) , one can also get the gauge coupling renormalization constant 
Z as from quark-gluino-squark vertex 

Z Qa = ZqgqZq %g Zq 1 (51) 

which, due to supersymmetry, of course coincides with the results of eqs. ff2"T]) and (jSjj). 



3.4 The e-scalar sector renormalization 

The renormalization of the ^-scalar sector in SUSY QCD goes along the same lines as for the ordinary 
QCD in dimensional reduction. Here we will remind the tricks used to obtain the corresponding 
renormalization constants and write down their values in SUSY QCD. The part of the SUSY QCD 
Lagrangian containing e-scalar fields could be obtained through the dimensional reduction of the 
Lagrangian from eq. (jlOj) . As a result we have 

£e = ^ (D^G^ D^Gcr — gz<ilaT a G a a q — ^glf abc f ade G b a G c a ,G^G e a , 

-i^T^g + glqT a T b G a a G b a q. (52) 

Here G a t a ') denote e-scalar fields, and indices a, a' belong to 2e subspace. 

A standard way is to use the Feynman rules derived from this Lagrangian to account for the 
contribution of e-scalars. However, it turns out that for some problems one can use far simpler 
arguments to get a desired solution. For example, we are interested in both full and separate 
contributions of vectors and e-scalars to the quantity, like a Green function or amplitude, without 
external vector indices. So, vector and e-scalar fields are located on internal lines only. Then to 
solve this problem one should 

(a) perform (i-dimensional vector and spinor algebras 

(b) for full contribution - replace d in the numerator of scalarized expression with 4; 
for e-scalar contribution - replace d n with 4 n — (4 — 2e) n . 
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q A\ q 



q ^ y g ^ y e ~\ y q 
Figure 5: One- loop diagrams for e-scalar Yukawa coupling renormalization constant. 





When we deal with a Green function containing external vector indices and want to extract a 
contribution of either external vector or e-scalar fields, the corresponding projector onto 4 — 2e or 
2e subspaces should be constructed. Internal e-scalars and vectors are treated in the same way, as 
described above. At this point one should keep in mind that particle momenta are always orthogonal 
to 2e subspace. For example, in the case of the e-scalar propagator (see FigHJ) we use ^g^ u (the 
e-scalar propagator is always diagonal, see the notice above) as a projector and after contraction 
of external indices replace d n with 4 n — (4 — 2e) n . At the same moment, evaluating corrections 
to the interaction between the e-scalar and quarks (see FigEJ), required for the determination of 
the Yukawa charge renormalization, one should follow the same steps, as in the case of e-scalar 
propagator with the only difference that now the projector is given by 



Tr 



1 



8C fCa£ 



fl } 5 



(53) 



where V 1 is an interaction vertex of the gluon field with quarks computed in SUSY QCD without 
separating e-scalar contributions, but using the 4-dimensional vector and spinor algebras. 
In the one-loop order the e-scalar field renormalization constant is given by 



1-^(6-^(2-0)-- 

47T v ' 6 



(54) 



For the one-loop Yukawa charge renormalization constant, describing the interaction of the e-scalar 
with fermions in SUSY QCD, we have 



'Y 



^qsq^q Z s 



Y- (3^4 - 6) -. 
47r e 



(55) 



Note that in SUSY QCD, contrary to the case of nonsupersymmetric QCD, it coincides with the 
gauge charge renormalization constant. It is due to supersymmetry, which now protects a tree- level 
coincidence of the £-scalar and vector coupling constants. 

The one-loop nonminimal mass counterterm for e-scalars reads 



5m„ = — 




+ 2 



-E 



m , 



+ 1 



+ Cj 



^ 2 



m~ 



+ 2 



(56) 



It is necessary to insure that the pole mass of the e-scalar is zero. Here rif denotes the sum over 
different quark flavours and nj is used to denote the sum over different squarks. It is just this 
counterterm that cancels m^ ard terms occurring in the two-loop diagrams and ensures decoupling 
of large scale physics for a physical quantity like a quark pole mass. A way of renormalizing the e- 
scalar mass nonminimally was first proposed in [7] and was shown to be a consistent renormalization 
procedure. In fact, there are two theoretically admissible prescriptions to deal with the e-scalar mass 
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renormalization. One is to renormalize the e-scalar mass minimally j2Hj and the other is the total 
subtraction of loop corrections to the e-scalar mass [7j. In the first approach, the e-scalar mass 
becomes a running quantity and physical quantities, like pole mass in our case, will in general 
depend on it together with other physical renormalized parameters. However, the e-scalar mass is 
not a physical quantity and it is more correct to treat it as an artifact of our regularization and 
renormalization schemes. From this point of view the second approach is more natural, if we are 
going to follow as close as possible the original idea of dimensional reduction suggesting the zero tree- 
level mass for both £-scalars and gauge bosons. Of course, all renormalization schemes are equivalent 
and could be related via a final recalculation of parameters. Introduction of nonminimal subtraction 
for the e-scalar mass leads to the fact that the bare quark mass is not anymore a renormalization 
group invariant. The point is that the physical quark mass depends on two mutually correlated 
bare masses: quark mass itself and e-scalar mass. Both these quantities are jl dependent and only 
the physical pole mass is RG invariant. For precisely this reason the RG equations for the bare 
quark mass should be written with great care. 



4 Results of calculation 

Here we present the results of our calculation. Corrections to the pole masses are parameterized as 
follows: 

44* = l + (^, (57) 
m V m J 

where Am starts from the one-loop order. 

Through calculations we reproduced all known results about radiative corrections to the pole 

masses of the bottom and top quarks: one-loop MSSM [TJl EH EH EE1 HEj and two loop QCD [HJ Ej 

(in MS ) and 7\ (in DR ) corrections. As an additional check, the renormalization group analysis can 

be applied to confirm independently the 1/e and 1/e 2 terms. We get complete agreement between 

our pole terms and those from RG for the DR quark mass renormalization constant. And the last 

but not least: calculations were performed in the linear gauge with a free QCD gauge parameter. 

The correction to the pole mass ([57)1 however appears to be gauge independent. 



4.1 One- loop result 



The corrections to the bottom and top quark pole masses in the one-loop order have been well 
known for a long time for QCD as well as for MSSM. Here we reproduce them only up to terms of 
the order 0(m 2 oft /m1 ard ) in a large mass expansion procedure. For pure QCD we have 



Am; 



QCD 



a. 




and the squark-gluino MSSM contribution is given by 



/ Ami 



+ 



+ 




2 2 

m~ — m~ 

02 & 



(5? 



(59) 
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which coincide with the results of IT3*| ITl} ITK| IT?)] . To obtain similar expressions for the case 
of the top one should perform an obvious substitution b — > t. The results for squark-chargino and 
higgs-quarks loops also known, but we are interested only in 0(a s ) and 0(a 2 s ) corrections in this 
paper. 



4.2 Two-loop result 



The two-loop QCD results in the DR-scheme for the 6-quark (with four light quarks) 



Ami 



m b ) 
+C F 



QCD 



c h 



OLsY 

Air) 



623 



8C 2 + ln 2 



mi 



rut 



^ + 30 C2 - 48 ln(2) & + 12 Ca + | hi I ^ ] i -hr j ^ 



mi 
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9, 



In 



mi 



mi 



mi 



+C A 



1 0Q3 1 7Q 

l " > -8G + 24 ln(2)C 2 -6C 3 -— In 



24 



6 



Hi 2 

H In 2 

2 



^ 2 



mi 



mi 



+ 26 In ^ - 6 In" 



and for the t-quark (with five light quarks) 
-Am^ QCD ^ 2 ' 



to* 
+C F 

+ C A 



a 
4irJ 



V _ 



43 - 12(2 + 26 In 



59 



^ + 30 C 2 - 48 ln(2) Ca + 12 C 3 + ~ In ( ?f ] + ^ln 2 ( ^ 



-6 In 2 



1 0Q3 1 7Q 

^-8C 2 + 24 ln(2)C 2 -6C 3 -^ln 



'ml 
^ 2 



n , 2 
H In 2 

2 



coincide with the first terms of the expansion in m^jvrit from [7] 

We find very a simple answer to one particular limit. Namely, in the case, when mr 



(60) 



(61) 



m 



tR 



mi 



"h — mr — 

"L OR 

definition of a, 



m 



m susY = M, the expression looks like (see eq. (|T2|) for the 



MSSM 



m n 



+C F 
+C A 



Ait 

Ml) 



23 13, 
In 

24 6 

175 41 



47 



+ iln 2 

2 



+ 20 In 



+ 



M 



72 6 
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'Ml) 
IW) 
( M 2 



'Ml) 

- 3 In 



+ 6 In 



M 2 ' 
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In 



M 2 



In 



77^ 
^ 2 



'M 2 ' 



M 
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3 ~ 



2 In 



M 2 
^ 2 
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In 



M 2 
> 2 



+ 3 In 



J' 



M 



8 ^ /M 2N 



(62) 



m? but not its 



Note, that we put equal to each other the left and right squark masses m^ R 
physical mass eigenstates m? = m? . The latest is true only when a q = because a minimal mixing 



takes place when 6 q 



i~ = m 2 
<?i 

7r/4 and ra~ 



2 

91,2 



M 2 ± a g m g (for a 9 < 0). 
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Since the complete formulae 8 are too large to be presented here, we present only the influence of 
our results on the masses of b- and t-quarks and the spectra of supersymmetric particles obtained 
as solutions of RG equations. 

We incorporated our formulae into the SoftSUSY code |Hj and analyzed the difference in the 
predictions for particle masses obtained with and without our corrections to the b- and t-quark 
pole masses. To measure quantitatively influence on the spectra let us introduce for each species of 
particle the quantity 

^CODE _ m CODE+Corrections 
A P = m CODE ' ( 63 ) 

where p stands for quarks or different supersymmetric particles. Below we put in our analysis the 
renormalization scale p, to be equal to mz- 

In FiglHl we present the variation of A p from the above formula as a function of tan/3 for 
m = 400 GeV, m l/2 = 400 GeV, A = and M GUT = 1.9 x 10 16 GeV. We plotted A p for five 
different cases: 

• Bottom 1L — difference between the results obtained with code including and excluding 
MSSM one-loop corrections to the pole mass of the 6-quark 

• Top 1L — the same as above for the t-quark 

• Bottom — difference between the results obtained with the code containing two-loop SUSY 
QCD correction to the 6-quark together with MSSM one-loop corrections to the pole masses of b- 
and t-quarks and results obtained with the code containing MSSM one-loop corrections to the pole 
masses of b- and t- quarks 

• Top — the same as above for the t-quark 

• All — difference between the results obtained with the code containing two-loop SUSY QCD 
correction together with MSSM one-loop corrections to the pole masses of b- and t-quarks and the 
results obtained with the code without any MSSM (one-loop and two-loop) corrections to the pole 
masses of b- and t-quarks. 

From this figure one can see that two-loop correction to the relation between the 6-quark pole and 
DR masses gives a sizeable contribution only to the solution obtained with the SoftSUSY program 
for the DR 6-quark mass at mz and almost does not influence spectra of SUSY particles. However, 
if one considers large tan j3 it becomes important for other particle spectra and thus should be taken 
into account in an accurate RG analysis. Note that our two-loop SUSY QCD correction for a wide 
range of parameter space is always less than a similar one-loop contribution. The two-loop MSSM 
correction to the t-quark pole mass, on the other hand, contributes to almost all SUSY particle 
masses obtained as a solution of renormalization group equations with universal boundary conditions 
at the GUT scale. This 2-loop correction is smaller than 1-loop MSSM contribution, however its 
effect on particle spectra is greater then corresponding 1-loop and 2-loop MSSM contributions to 
the 6-quark pole mass almost everywhere in MSSM parameter space. This could be attributed to 
potentially large corrections coming from higher terms in the mt/ Msusy expansion. 

To see the dependence of A p on mo, mi/2 and A , we use so called a set of benchmark points 
and parameter lines in MSSM parameter space from [23 123 EE I2H] which corresponds to different 
scenarios in the search for Supersymmetry at present and future colliders. While a detailed scan 
over the more-than-hundred-dimensional parameter space of MSSM is clearly not practicable, even a 
sampling of the three- (four-) dimensional parameter space of uiq, and A (tan/?) is beyond the 
present capabilities of phenomenological studies, especially when one tries to simulate experimental 
signatures of supersymmetric particles within detectors. For this reason, one often resorts to specific 



3 They could be found in the source archive file for this paper submitted to the HEP database 
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benchmark scenarios, i.e., one studies only specific parameter points or at best samples the one- 
dimensional parameter space (the latter is sometimes called a model line |28j), which exhibit specific 
characteristics of MSSM parameter space. 

Some recent proposals for SUSY benchmark scenarios may be found in (2E1 123 EH]- We refer 
to the "Snowmass Points and Slopes" (SPS) [231211 which consists of model lines ("slopes"), i.e., 
continuous sets of parameters depending on the one dimensionful parameter and specific benchmark 
points, where each model line goes through one of the benchmark points. 

In Figd we present the results for A p as functions defined along Model Line A 

m = —A = 0.4mi/2, 77H/2 varies, tan/3 = 10, \i > 0. 

where the benchmark point is 

m = 100 GeV, m 1/2 = 250 GeV, A = -100 GeV, tan /3 = 10, p. > 0. 

Here and below /1 is the supersymmetric Higgs mass parameter 

As one can easily see from FigJTJ varying mi/2 at fixed tan/? 2-loop SUSY QCD correction to 
the 6-quark pole mass almost does not contribute to heavy supersymmetric particle spectra. On 
the other hand, our correction to the i-quark pole mass is comparable to the corresponding 1-loop 
contribution. 

In Fig|Hl we present the results for A p defined as functions along Model Line C 
m = mi/2, mi/ 2 varies, A = 0, tan/3 = 35, \i > 0. 
where the benchmark point is 

m = 300 GeV, m 1/2 = 300 GeV, A = GeV, tan/? = 35, /i > 0. 

Situation with Fig|S]is analogous to that of FigJTl even if we increase tan/? from 10 up to 35. 
Importance of our corrections increases only at very large values of tan/3 (50 and more), as one can 
see from FigEl • 

For such a large tan/3 we just used some particular values m and from the allowed region 
from [3U]. We use the following parameterization: 

m = mi/2, mi/2 varies, A = 0, tan/3 = 55, p > 0. 

The results for A p are shown in FigHJ and we can see that in this case our 2-loop correction to the 
6-quark pole mass gives a sizeable contribution to the masses of gluino and sbottoms. 

Now let us discuss our limiting formula eq. ([62)1 . We have studied the difference between our 
full result and this approximate expression in different regions of MSSM parameter space using the 
same strategy as above. As a result we find (see Fig. [TUJ) that for the 6-quark this difference does 
not exceed 12% for the parameters region considered above, both for the dependence on tan (3 and 
77ii/2 . For the t-quark this is even better. 
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5 Conclusion 



We presented the results for two-loop MSSM corrections to the relation between pole and running 
masses of heavy quarks up to the 0(a 2 s ) order. We provided a detail analysis of the value of these 
corrections in different regions of parameter space and discussed their impact on SUSY particle 
spectra. We would like to note that our results presented here may also be included in the codes 
of programs like Isajet [3T], SuSpect [32J and used for predictive phenomenological analysis like 
[3*0"] . [""3~| . "3*"Ij . jUSj. In one of our next papers we suppose to calculate missing corrections from 
stop-chargino loops to the 6-quark pole mass and to provide more terms in the expansion in the 
relation between the top pole and DR masses. 

The authors would like to thank K.Chetyrkin, M.Kalmykov, D.Kazakov and V.Smirnov for 
fruitful discussions and multiple comments. This work of A.O. and O.V. was supported by DFG- 
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264/2-1) and by BMBF under grant No 05HT9VKB0. Financial support for A.B. and V.V. from 
RFBR grants # 02-02-16889 and # 00-15-96691 is kindly acknowledged. 
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Figure 6: A p from the formula (J63j) as functions of tan (3. In the first picture we plotted a spectrum 
of masses (see also explanation in the text). 
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Figure 7: A p from the formula (|63j) as functions defined along Model Line A from corrections 
depending on mi/ 2 . In the first picture we plotted a spectrum of masses 
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Figure 8: A p from the formula (J63)) as functions defined along Model Line C from corrections 
depending on ra\j%. In the first picture we plotted a spectrum of masses 
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Figure 10: Difference between A p with the full two- loop results and the limit from eq.([62|). Upper 
plots corresponds to the t-quark and bottom ones are for the 6-quark. Left column shows the 
dependence on m 1( / 2 of the difference between A p and A^ 111 ^. Right column shows the relative 

dependence on tan (3 of the difference between A p and A^ m ^ for the same parameters as in text. 
Line A and Line C correspond to Model Line A (see Fig. E|) and Model Line C (see Fig. |HJ) 
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